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Abstract. This paper addresses the synchronized region problem, which is reduced to a matrix 
stability problem, for complex dynamical networks. For any natural number n, the existence of a 
network which has n disconnected synchronized regions is theoretically demonstrated. This shows 
the complexity in network synchronization. Convexity characteristic of stability for matrix pencils 
is further discussed. Smooth and generalized smooth Chua's circuit networks are finally discussed 
as examples for illustration. 
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1 Introduction and problem formulation 

The subject of network synchronization has recently attracted increasing attention from various 
fields (see [H El El [H El HH H21 HH [201 E21 [231 [2H EH1 [26] and references therein). Of particular 
interest is how the synchronization ability depends on various structural parameters of the network, 
such as average distance, clustering coefficient, coupling strength, degree distribution and weight 
distribution. Some important results have been established for such problems by introducing the 
notions of master stability function and synchronized region [H [HI [T2"l [13"1 [TEl 127j . It is natural 
to expect strong synchronization ability at small cost [13] • In fact, a key factor influencing the 
synchronization ability is the characterization of the network synchronized region, as studied in 
[HI 1101 [To]. Obviously, the larger the synchronized region, the easier the synchronization. Some 
examples for the existence of two and three disconnected synchronized regions are demonstrated in 
|10j . This paper attempts to explore the existence of multiple disconnected synchronized regions 
for various complex dynamical networks. 
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Consider a dynamical network consisting of N coupled identical nodes, with each node being an 
n-dimensional dynamical system, described by 

N 

&i = f(xi) + c'^2a ij H(x j ), i = 1,2, ■■■ ,N, (1) 
3=1 

where Xi = (xn, Xii, ■ ■ ■ ,Xin) S K n is the state vector of node i, /(•) : M. n — » 1" is a smooth 
vector-valued function, constant c > represents the coupling strength, H(-) : R n — > R n is called 
the inner linking function, and A = (aij)NxN is called the outer coupling matrix, which represents 
the coupling configuration of the entire network. Generally, A is an irreducible matrix, and if the 
entries of A satisfy 

N 

<Hi = — ^ aij, i = 1,2,- ■■ ,N, 

then network (JTJ) is called a diffusively coupled network. In this case, zero is an eigenvalue of A 
with multiplicity 1 and all the other eigenvalues of A are strictly negative, which are denoted by 

= Ai > A 2 > A 3 > • • • > X N . (2) 

The dynamical network ([T]) is said to achieve (asymptotical) synchronization if 

xi(t) — > X2{t) —>•••—> X]y(t) — > s(t), as t — > oo, (3) 

where, because of the diffusive coupling configuration, the synchronous state s(t) G M n is a solution 
of an individual node, i.e., s(t) = f(s(t)). Here, s(t) can be an equilibrium point, a periodic orbit, 
or even a chaotic orbit. 

As shown in [Eld!], the stability of the synchronized solution xi(t) = x 2 (i) = • • • = x^(t) = s(t) 
can be determined by analyzing the following equation, known as the master stability equation: 

u = [Df(s(t)) + aDH(s(t))]LJ, (4) 

where a€l, and Df(s(t)) and DH(s(t)) are the Jacobian matrices of functions / and H at s{t), 
respectively. 

The largest Lyapunov exponent L max of network (pQ), which can be calculated from system @ 
and is a function of a, is referred to as the master stability function. In addition, the region S of 
negative real a where L max is also negative is called the synchronized region of network ([1]). Based 
on the results of |12[ Ho] , the synchronized solution of network ([T]) is asymptotically stable if, and 
only if, 

cA fc ES, fc = 2,3,--- ,N. (5) 

The synchronized region S can be an unbounded region, a bounded region, an empty set, or 
a union of several regions. If the synchronous state is an equilibrium point, then Df{s{t)) and 
DH(s(t)) reduce to constant matrices, denoted by F and H, respectively. In this case, system (J3J) 
becomes 

uj = [F + aH]uj. (6) 



Hence, the synchronized region S becomes the stability region of F + aH with respect to parameter 
a. This paper mainly studies this case when the synchronous state is an equilibrium point. 

The rest of this paper is organized as follows. In Section 2, the disconnected stability region 
problem for the matrix pencil F + aH is studied, where the existence of multiple disconnected 
stability regions is theoretically proved. In Section 3, the characteristics of matrix convexity for 
the stability of matrix pencils are discussed, where some conditions for testing the stability or 
instability of convex combinations of two vertex matrices are established. In Section 4, smooth 
Chua's circuit networks are simulated to illustrate the theoretical results. The paper is concluded 
by the last section. 



2 Disconnected stability regions for matrix pencils 

As discussed in the previous section, when the synchronization state is an equilibrium state, the 
synchronized region problem reduces to a stability problem of the matrix pencil F + aH with 
respect to parameter a. In this section, the characteristics of disconnected stability regions for 
such matrices are studied. In order to discuss this problem in the real parameter domain, the 
following lemmas are necessary. 
Lemma 1 If the real polynomial 

p(s) = s n + Tn-is™" 1 + • • • + 7i s + 70 (7o > 0) 

is stable, then for any scalar e, < e < 70, the following polynomial 

Pe (s) = s n + 7n-is n_1 + • • • + lis + e 

is stable. 

Proof Given that p(s) is stable, polynomial p e (s) is stable if, and only if, p(s) — e is stable for 
all < e < 70, or equivalently, the function 

e 

p(s) 



p(sj 

is stable. Further, this is equivalent to that the Nyquist plot of — does not enclose the point 
(—1,0) for all < e < 70, which obviously holds. □ 

Lemma 2 Given a polynomial p(a) = (a + l)(a + 2) • • • (a + n) with variable a and n > 2, 
there is a scalar (3 > such that p{ct) — (3 n has n negative real roots. 

Proof Take (3 > such that (3 n < ±(0.5 x 1.5 x • • • x ([§] - 0.5)). Then, one can get 

p(0)-f3 n > 0, p(-l)-P n < 0, p(-2.5)-p n > 0, • ■ ■ , p (-2 x [|] - 0.5) -f3 n > 0, p(-n)-[3 n < 0. 

Therefore, the sign of p{a)—(3 n changes n times on the negative real axis. This means that p(a) — (3 n 
has n real roots on the negative real axis. □ 
Lemma 3 Given two scalars (3q and (3 with (3 > and (3 — (3q > 0, there are scalars < a\ < 
■ ■ ■ < a n such that «iQ2 • • • a n = (3 and all roots of p(a) = (a + a\)(a + 02) ■ • • (a + oi n ) — {(3 — (3q) 
are real. 



Proof By the method of Lemma 2, it suffices to prove this lemma by choosing a% with the above 
constraints such that the sign of p(a) changes n times on the real axis. □ 

With the above lemmas, one can get the following results. 

Theorem 1 For any natural number n, there are matrices F and H of order 2(n — 1) such that 
F + aH has n disconnected stable regions with respect to parameter a. 
Proof As shown in Lemma 2, one may take (3 > such that 



p(a) = (a + l)(a + 2) • • • (a + 2(n - 1)) - 0^ n ~^ = 
has 2(n — 1) real roots, denoted by /3i, $2, ■ • ■ , fh(n-l)- Then, take 

/ -1 ••• \ / 0i ■■■ \ 



(7) 



H 




\1 



-1 

J 



F, 



'•• /3 2 n-3 

V -/9 2 (n-l) • • • / 



and F = — Pl^in-l) + w here J2(n-l) ^ s the identity matrix of order 2(n — 1). Obviously, the 
characteristic polynomial of F + CK-ff is 

det(sJ -F-aH) = (s + /J) 2 ^ 1 ) + (a - /?i)(a - /3 2 ) • • • (a - /3 2(n _i)). 

Using ((TJ), one has 

det(sl — F — aH) = (s + p) 2 ^-^ - p 2(n -^ + (a + 1)(q + 2) ■ • • (a + 2(n - 1)). 

The constant term in det(s/ — F — aH) is 

(a + l)(a + 2)---(a + 2(n- 1)), 

which is larger than zero if the parameter a is located in the following n regions: 

(0, -1), (-2, -3), • • • , (-2(n - 2), -2n + 3), (-2(n - 1), -00), (8) 

and is smaller than zero if a is located in the following n — 1 regions: 

(-1, -2), (-3, -4), • ■ ■ , (-2n + 3, -2(n - 1)). 

Obviously, by Lemma 1 det (si — F — aH) has n disconnected stable regions with respect to 
parameter a, which are contained in the n regions shown in 0, respectively. □ 

Combining with the discussions in Section 1, for any natural number n, Theorem 1 shows the 
existence of a network which has n disconnected synchronized regions. However, for a general 
network, the node equation is given, i.e., F is given, which can not be chosen arbitrarily. In this 
case, one may apply the following result with a chosen inner linking matrix H. 

Theorem 2 For any given real stable matrix F of order n, suppose det(sI-F) = s ra +7 n _is n_1 + 
• • • + 71 s + 70, and every eigenvalue of F corresponds to only one Jordan form. If there is a scalar 
Po ^ such that p(s) = s n + 7„_is n ~ 1 + • • • + 71s + 70 — fto is stable and p(s) has ni pairs of 



conjugate complex eigenvalues, then there exists a real matrix H such that F + aH has [ n 2 "' ] + 1 
disconnected stable regions with respect to parameter a. 
Proof First, suppose that there is a scalar 0q such that 

p(s) = s n + 7n-is n_1 + • • • + 71 s + 70 - A) 

is stable, with n real roots denoted by Aoi, • • • , Xon- Following Lemma 3, take scalars < a\ < 
■ ■ ■ < a n such that aia 2 • • • a n = 70 and all roots of (a + ai)(a + a 2 ) • • • (a + a n ) — (70 — /?o) are 
real, denoted by — 0i, • • • , — n . Consequently, 

p(a) = (a + ai)(a + a 2 ) • • • (a + a n ) - (70 - A)) = (a + /?i)(a + • • • (a + /?„). (9) 

Obviously, 0i ■ • ■ n = 0o- Furthermore, take 

/ 1 • • • \ / Aoi 0\ ■ ■ ■ \ 



F = 

! ■■■ /? n _i 

/ V -/?„ • • • A n / 

(s-A „)+A) = s n +7 n _is n " 1 + - • -+7is+7o = det(sI--F). Hence 





V-i 

Then, det(sI-F ) = (s-A i) 
Fq is similar to F, since each eigenvalue of F corresponds to only one Jordan form. Moreover, 

det(s/ - F - aH ) = (s - A i) • • • (s - A „) + (a + (5i) ■ ■ ■ (a + /3 n ) 

= (s - Aoi) •••(«- A 0n ) - (70 - A)) + (a + ai)(a + a 2 ) • • • (a + a n ). 
The constant term in det(s/ — Fq — aH®) is 

(a + qi)(q + a 2 ) • • • (a + a n ), 
which is larger than zero if the parameter a is located in the following [5] + 1 regions: 

(0, -qi), (-a 2 , -03), • • • , (10) 
and is smaller than zero if a is located in the following n — [^] regions: 

(-ai,-a 2 ), (-a 3 ,-a 4 ), 

Thus, by Lemma 1, Fq + ckHq has [^] + 1 disconnected stable regions, which are located in the 
regions shown in (|10p . Since -Fo is similar to F, there exists a nonsingular matrix P such that 
P^ 1 FqP = F. Therefore, H = P^ 1 HqP is the matrix to be found. And F + aH has the same 
stable regions as Fo + oHq. 

Then, assume that there are some conjugate complex pairs in Aoi,-- - ,Ao n - For simplicity, 
suppose that there is only one pair of conjugate complex eigenvalues, Aoi = £1 + Vih ^02 = £1 — Vih 
and A03, • • • , Ao n are all real. 

Similarly to the above proof, take scalars < a 2 < ■ ■ ■ < a n such that 02^2 • • • & n = 70 and all 
roots of (a + a 2 )(a + 03) •••(« + a n ) — (70 — po) are real, denoted by — /3 2 , • • • , —(3 n . Consequently, 



p(a) = (a + a 2 )(a + a 3 ) • • • (a + a n ) - (70 - 0o) = (a + ^2) (a + 02,) • • • (a + ^n)- (H) 



Obviously, ft ■ ■ ■ ft = ft- Furthermore, take 
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Obviously, det(sI-Fo) = (s-2£is+ff+7/f )(s-A 03 ) ■ ■ ■ (s-A n)+ft = s n +7n-is n_1 +- ■ -+7is+7o = 
det(s/ — F). Hence, Fo is similar to F. Moreover, 

det(sl -F - aH ) = (s - 2^s + g + rfi )(s - A 03 ) ■ • • (s - A 0n ) + (a + ft) • • • (a + ft) 

= s n + Tn-is™- 1 + . . . + 7lS + ( a + a 2 ) ■ ■ ■ (a + a„). 
The constant term in det(s/ — Fq — aH$) is 

(a + a 2 ) • • • (a + ««), 



which is larger than zero if the parameter a is located in the following [^^rp] + 1 regions: 

(0, -a 2 ), (-a 3 , -a 4 ), 



(12) 



Repeating the process as above, one can complete the proof easilly. □ 

Remark 1 For simplicity, in both Theorems 1 and 2, the parameter a only appears in the 
constant term of the characteristic polynomial of F + aH. In this way, the constant term in 
det(sl — F — aH), which is a polynomial of parameter a, simply determines the disconnected 
stable regions. If a appears in the higher-order terms of det(sl — F — aH), the problem becomes 
harder to solve, leaving an interesting topic for further research. 

Remark 2 In order to guarantee H be a real matrix, two cases are considered in the proof of 
Theorem 2, i.e., there are or there are no conjugate complex pairs in Aoj, i = 1, • • • , n. If H can be 
chosen to be a complex matrix, the proof of Theorem 2 will be simplified and H may be chosen 
such that F + aH has [^] + 1 disconnected stable regions. In addition, if all \oi,i = 1, • • • , n, are 
complex scalars, then there exists a real H such that F + aH has at least [j] + 1 disconnected 
stable regions. 

According to the above discussions, one can also choose a suitable H such that F + aH has only 
one convex stable region with respect to parameter a, as further discussed below. 



3 Characteristics of convexity for stability of matrix pencils 

In the previous section, it shows the existence of any n disconnected stable regions of the matrix 
pencil F + aH. Contrary to this non-convexity, given two parameter values a\ and q 2 , whether 
or not the stability of F + a\H and F + a 2 ^T implies the stability of F + {\a\ + (1 — A)a: 2 )il, 



for all < A < 1, is an interesting problem. Obviously, a good understanding of this convexity 
characteristic is useful for enhancing the stability of the matrix pencil F + aH . 

Lemma 4 Suppose that F + a\H and F + a 2 H are stable, and the rank of H is 1. Let 
H = be, where b is a column vector and c is a row vector with compatible dimensions, and (F, b) 
be controllable. Then the following conditions are equivalent to each other: 

(i) \{F + aiH)- 1 + (1 - \)(F + a 2 H) is stable for all < A < 1. 

(ii) There is a common matrix P = P T such that 

P{F + ai H) + (F + ai H) T P <0, i = 1,2. 

(iii) (F + a\H){F + a 2 H) does not have negative real eigenvalues. 

(iv) 1 - Re{(« 2 - a!)c(jwl -F- > 0, Vu) G R. 
Further, if any one of (i)-(iv) holds, one has 

(*) X(F + ai#) + (1 - X)(F + a 2 #) is stable for all < A < 1. 

Proof See |18l 119] for the equivalences among (ii)-(iv). Now, if (i) holds, then 

det(A(F + a^y 1 + (1 - \){F + a 2 H)) / 0, V < A < 1. 

Equivalently, 

det (j—J 1 + ( F + a i R )( F + a 2 H)J + 0, V < A < 1, 

which implies (iii). On the other hand, by a simple congruence transformation with matrix (F + 
ot\H)~ l } (ii) implies the existence of a common matrix P = P T such that 

P(F + a^y 1 + (F + ai H)- T P < 0, P(F + a 2 H) + (F + a 2 H) T P < 0, (13) 

which implies (i). And, obviously, (ii) implies (*). This completes the proof. □ 
Remark 3 Any one of Lemma 4 (i)-(iv) implies (*). However, generally, (*) does not imply 
the other conditions of Lemma 4. For example, with 



F = 







' 






,H= ( 




: 











, ai = 0, a 2 = 0.9, 



obviously Lemma 4 (*) holds for the above matrices, but all Lemma 4 (i)-(iv) do not hold. 

Remark 4 Obviously, Lemma 4 (ii) is equivalent to (|13|) . However, since the rank of (F + 
ai_£f) -1 — (F + a 2 H) is generally not 1, (fT3|) is not equivalent to that (F + ot\H)~ l {F + 02-?/) does 
not have negative real eigenvalues. In fact, (F + aiHy 1 (F + a 2 H) = (F + aiff) _1 (i ? + + 
(«2 — ol\)H) = I + («2 — ai)(i ? + a\H)^ 1 H. That (F + a\H)^ l {F + Q2-ff) does not have negative 
real eigenvalues is equivalent to 1 + (02 — ai)c(F + a\H)~ l b > 0, which is Lemma 4 (iv) at jw = 0. 
The example given in Remark 3 can also be an example for this. 

Remark 5 It is not necessary to require P > in Lemma 4 (ii). The positive definiteness of P 
is naturally guaranteed by the stability of F + aiH, i = 1,2. For this property, the following lemma 
is useful. 



Lemma 5 |9j Suppose that F and P = P T are matrices of order n. If PF + F T P < 0, then F 
has no eigenvalues on the imaginary axis, det(P) 7^ 0, and the number of positive eigenvalues of P 
is equal to the number of eigenvalues of F with negative real parts. 

Corollary 1 Suppose the rank of H is 1 and F + ot\H is stable. Let H = be and (F, b) be 
controllable. Then, F+aH is stable for all a 6 (—00, ai] if, and only if, He{c(jwl— F— aiH)^ 1 ^ < 
0, Vto G R. 

As to the instability of matrix pencils, the following result can be obtained from Lemmas 4 and 

5. 

Theorem 3 Suppose that F + ot\H and F + a 2 H are unstable and do not have imaginary 
eigenvalues, and the rank of H is 1. Let H = be and (F,b) be controllable. Then, the conditions 
of Lemma 4 (ii)-(iv) are equivalent to each other. Further, if any one of Lemma 4 (ii)-(iv) holds, 
then F + (1 - X)a\H + Xa 2 H is unstable for all A G (0, 1). 

Proof Since the Kalman-Yakubovich-Popov lemma also holds for unstable state matrices [17j . 
Theorem 5 in [6] yields the equivalence between (ii) and (iv) of Lemma 4. Since F + a\H and 
F + a 2 H do not have imaginary eigenvalues, Theorem 3.1 in |19] gives the equivalence between (ii) 
and (hi) of Lemma 4. 

In what follows, it is to prove that Lemma 4 (iv) implies the instability of F+a\H+\(a 2 — a\)H. 
Suppose that there is a Ao G (0, 1) such that F + a\H + Ao(a2 — a\)H is stable. Since matrix 
eigenvalues change continuously with matrix parameters, by the instability of F + a\H, there exists 
a Ai, < Ai < Ao such that F + a±H + Ai(a2 — a\)H has an imaginary eigenvalue jwo, i.e., 

det(jwoI — F — ot\H — Ai(«2 — ct\)H) = 0. 

Since F + ot\H does not have imaginary eigenvalues, the above condition is equivalent to 

det (J-I - (a 2 - ax)(jw I — F — aiH)- 1 ^ = 0, 

or 

(«2 — &i)c(jwoI — F — aiH)~ 1 b = 0, 

Ai 

which is contrary to Lemma 4 (iv). This completes the proof. □ 
Remark 6 The common Lyapunov matrix problem for stable matrix pencils was studied in 
\18\ 119]. Theorem 3 above generalizes the similar results to unstable matrix pencils. In fact, any 
one of Lemma 4 (ii)-(vi) guarantees that transferring an eigenvalue between the left-half and right- 
half complex planes is impossible. Therefore, when any one of Lemma 4 (ii)-(vi) holds, F + a\H 
and F + a 2 H must have the same number of eigenvalues with positive real parts. In addition, 
the stability of the convex combinations of (F + aiH)^ 1 and F + a 2 H is equivalent to any one of 
Lemma 4 (ii)-(iv). However, for instability, this equivalence does not hold; that is, the instability 
of (F + ctiH)~ l + \(a 2 — ai)H, A G [0, 1], does not necessarily imply any one of Lemma 4 (ii)-(vi). 
For example, with 



obviously, (F + a\H) 1 + A(a 2 — a\)H is unstable for all A G [0, 1], but any one of Lemma 4 (ii)-(iv) 
does not hold. 

Although the characteristics of convexity for stability or instability of the matrix pencil F + aH 
have been discussed when the rank of H is 1, it is still hard to decide the stability or instability 
of the convex combinations of F + a±H and F + a 2 -ff for a general H. The results of Lemma 4 
and Theorem 3 are directly related to the existence of a common matrix P for two vertex matrices. 
For a general H, this common-matrix method is very conservative. In this case, nevertheless, the 
following lemma provides a less conservative criterion |15| . 

Lemma 6 Suppose that F + a\H and F + a 2 ff are stable. If there are matrices Pi = Pf, P2 = 
P-f, G and V such that 



-G — G T Pi-V T + G(F + ai H) 

P l -V + {F + aiH) T G T V(F + c^H) + (F + ai.H) T V T 



< 0, i = 1,2, 



then F + Xa x H + (1 - X)a 2 H is stable for all < A < 1. □ 

In the above lemma, by introducing new slack matrices G and V, the symmetrical matrices Pi 
and P2 can be chosen parameter-dependent for the study of stability of the convex combination of 
F + a\H and F + a 2 H. Similarly to the method used in [3 [10] , one can also design controllers to 
enlarge stability regions by the above-discussed method. 

For instability of matrix pencils, one can likewise obtain the following result. 

Theorem 4 Suppose that F + a\H and F + a2H are unstable. If there are matrices P\ = 
Pf, P 2 = P 2 T , G and V such that 



—G — G T Pi-V T + G(F + ai H) 

Pi-V + {F + ai H) T G T V(F + ai H) + (F + ai H) T V T 



< 0, i = 1,2, 



then F + (A«i + (1 - \)a 2 )H is unstable for all < A < 1. 

Proof For any A G (0, 1), taking a convex combination between two inequalities in Theorem 4, 
one gets 

( — G — G T XP 1 + (1-X)P 2 -VT + GF X \ 

{XP 1 + (1-\)P 2 -V + FTG T VF x + FlV T ) ' [ ' 

where F\ = F + (Acti + (1 — \)a2)H. Then, by the similar method used in [T2] , the above inequality 
is equivalent to 

(APi + (1 - A)P 2 )P A + P A T (APi + (1 - A)P 2 ) < 0. 

Here, F\ must be unstable. If it was stable, then, as proved in Theorem 3, there would exist 
Ao G (0, A) such that F\ = F + (Ao«i + (1 — Xo)a 2 )H has an imaginary eigenvalue jwq. Then, by 
Lemma 5, 

(A Pi + (1 - A )P 2 )Pa + ^a (A Pi + (1 - A )P 2 ) 
can not be strictly negative definite, which is contrary to (|14|) . □ 

Remark 7 Theorem 4 generalizes the method of [15J to the instability of matrix pencils. Obvi- 
ously, the instability of the matrix pencil F + aH is important in desynchronization problems. As 



discussed in Remark 6, if Theorem 4 holds, transferring any eigenvalue of F + a\H and F + a 2 H 
between the left-half and right-half complex planes is impossible. Therefore, F + a±H and F + a 2 H 
must have the same number of eigenvalues with positive real parts when Theorem 4 holds. 



4 Synchronization of smooth Chua's circuit networks 



The vector X{ in ([I]) is (xn,Xi 2> Xis) 



Fx,. 



(16) 



In this section, consider the synchronization problem in a network of smooth Chua's circuits. 

Example 1 Consider the network (TjQ) consisting of the third-order smooth Chua's circuits 
in which each node equation is 

xn = —kaxn + kaxi2 — ka(axf x + bxn), 

x i2 = kxa - kx i2 + kx i3 , (15) 
£i3 = -k(3x i2 - k^x i3 . 

■ T here. Linearizing (|15p at its zero equilibrium gives 

— ka — kab ka \ 
F= k -k k 

I —kj3 -kj J 

Take k = l,a = —0.1,/? = —1,7 = l,a = 1,6 = —25. Then F is stable, i.e., the node system 
(I15p is locally stable about zero. One can easily take a parameter (3q = —0.8 such that all roots of 
det(sl — F) — (3q are real. Following the method of Theorem 2, take a\ = 0.01, a 2 = 1,03 = 10, 
and 

0.8348 9.6619 2.6591 
0.1002 0.0694 0.1005 
-0.3254 -8.5837 -0.9042 
Then, by simply computation, one knows that F + aH has two disconnected stable regions: S\ = 
[—0.0099,0] and S2 = [—2.225,-1). Therefore, the entire synchronized region is SiljS^. Further, 
suppose that the number of nodes is N = 10, and the outer coupled matrix A is a globally coupled 
matrix, i.e., all the diagonal entries of A are —9 and the other entries are 1, which has eigenvalues 

A 1 = 0,A 2 = --- = A 10 = -10. 

Then, by ([5|), network (TjQ) with the above parameter values achieves local synchronization when 
the coupling strength c satisfies c G [0,0.00099] or c S (0.1,0.2225]. Figures 1 and 2 show the 
synchronization and non-synchronization behaviors of this network. 
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(a) c = 0.0005 6 [0, 0.00099]. (b) c = 0.2 G (0.1, 0.2225]. 

Fig. 1 Synchronization of network ([1]) with different coupling strengths. 




(a) c = 0.02 G (0.001, 0.1). (b) c = 0.3 G (0.2225, +oo). 

Fig. 2 Non-synchronization of network (pQ) with different coupling strengths. 

Example 2 Consider the network (p} consisting of the forth-order generalized smooth Chua's 
circuits [7J, with node equation 



in = —kaxn + kax i2 — ka(axf 1 + bxn), 
x i2 = kxn - kx i2 + kx i3 , 
x i3 = kf3x i2 + kjx i4 , 

Xi4 = -0.1x i2 . 

The vector Xi in JT} is (xn, Xi2, Xis, Xi<±) T here. Linearizing p!7|) at its zero equilibrium yields 



(17) 



Fx i} F 



( 



\ 



-ka — kab 
k 





ka 
-k 

k/3 





k 



\ 




kj 



-0.1 



(18) 



/ 



Take parameters k = 3, a - 0.1, ^ = -0.2,7 = 0.2, a = 1 and b = -25. Then F in CEHD is 
stable, i.e., the node system (|17p is locally stable about zero. One can easily take a parameter 
(3q = —0.1 such that all roots of det(s/ — F) — (3q are real. Following the method of Theorem 2, 
take a\ = 0.1, a 2 = 0.5, «3 = 2, 04 = 12.96, and 



H 



V 



0.8442 
-12.7738 
-10.3570 
-4.7028 



0.6319 0.3547 

-9.9676 19.4669 

-8.3421 18.4474 

-3.8156 8.5403 



-1.8905 \ 
-20.2986 
-20.5913 
-9.3240 J 

Then, by simply computation, one knows that F + aH has three disconnected stable regions: 
Si = (—0.1,0], S2 = (—2,-0.5) and S3 = [—12.95,-12.94], so the whole synchronized region is 
Si U 5*2 U S3. Theorem 2 implies that, the region £3 should be contained in (—00, —12.96), but due 
to the computing error S3 becomes [—12.95, —12.94], slightly off-set from the theoretical prediction. 

Similarly to Example 1, if N = 10 and the outer coupled matrix A is a globally coupled matrix, 
then both synchronization and non-synchronization phenomena can be discussed. 



5 Conclusion 



In this paper, the problem of disconnected synchronized regions has been carefully studied. When 
the synchronization state is an equilibrium point, the problem is reduced to the stability problem of 
matrix pencils. The existence of multiple disconnected synchronized regions is theoretically proved 
for network with higher-dimensional nodes. Further, the characteristics of convexity for matrix 
pencils has been discussed. Some test conditions for stability and instability of convex combinations 
of vertex matrices have also been established. Finally, networks of smooth and generalized smooth 
Chua's circuits have been simulated to illustrate the analytic results. 
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